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Abstract. Measurements of dijet cross sections for virtual photons are presented as a 
function of x® BS , the fraction of the virtual photon energy invested in the production 
of the dijet system, using the ZEUS detector. Comparisons to QCD predictions show 
that a resolved photon component is needed to describe the data up to values of the 
photon virtuality comparable to the scale of the interaction. 



INTRODUCTION 

It has been long established that real photons have a partonic structure from 
measurements of the photon structure function FJ in cy interactions [1] and ob- 
servation of resolved photon processes in 7p interactions [2] and of single- and 
double-resolved processes in 77 interactions [3]. Thus, it is natural to expect that 
virtual photons also have a partonic structure [4-7]. QCD predicts that the parton 
densities of virtual photons become logarithmically suppressed as the virtuality of 
the probed photon P 2 increases for fixed /1 2 , the scale of the interaction; fi 2 is 
usually taken as Q 2 (Q 2 is the virtuality of the probing photon) in deep inelastic 
scattering (DIS) &y and the jet transverse energy (E^) in jet production. 

The virtual photon structure function F^ at leading order (LO) in QCD is given 

by 



Fi (x 7 *,P 2 ,/i^J — ^^x 7 * e 2 [fq/Y (x 7 *, P 2 , ) + /g/ 7 *(x 7 *, P 2 , fi F ^ )], 

1 

where x 7 * is the fraction of the virtual photon momentum taken by the interacting 
parton, \j? F t is the virtual photon fragmentation scale and e q is the quark charge. 
The virtual photon quark densities / g / 7 * contain two terms, as in the case of the 
real photon, 
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FIGURE 1. (a) Resolved and (b) direct virtual-photon processes. 



fi/r* {x 7 * , P 2 , ^f 7 „ ) — fi/y- ( X T > P' \ ) + fi/r* m ( x Y i P 2 -, A^* )> 

one term associated to the non-perturbative hadronic component (/ had ) and a term 
/ anom , unique to the photon, which expresses the direct coupling of the photon 
to a qq pair, calculable in perturbative QCD (pQCD). The hadronic component 
is the one expected to decrease as P 2 increases. The first measurements of the 
structure function of the virtual photon were done by PLUTO [8]. Only recently 
new measurements from LEP have become available [9]. 

At HERA, the virtual photon structure is studied in jet production mediated by 
virtual photons. In LO QCD, two processes are expected to contribute to the jet 
cross section: resolved processes (figure la) in which the virtual photon interacts 
with the proton via its hadronic component giving two jets in the final state, and 
the direct processes (figure lb) in which the virtual photon interacts as a point-like 
particle with a parton from the proton. 

The dijet production cross sections in LO QCD for direct and resolved processes 
are given by 

a L °' ep = J dtt f r/e (y,P 2 ) fj/p(x p ,fi 2 Fp ) da{-f*j -> jet jet) 

and 

a LO,e P = J dn f r/e ( yi p*) / i/r (x 7 »,P 2 ,/4 7 J fi/ p {x p ,n 2 Fv ) da(ij jet jet), 

where the integrals are performed over the phase space represented by "dfi"; 
fr* /e(y> P 2 ) i s flux of virtual photons in the positron and y is the fraction of the 
positron energy taken by the virtual photon; fj/ p (x p , fip p ) are the parton densities 
in the proton and x p is the fraction of the proton momentum taken by parton j; 
lx 2 Fp is the proton fragmentation scale; and da(j*{i)j — > jet jet) is the subprocess 
cross section, calculable in pQCD. In the case of resolved processes, there is an ad- 
ditional ingredient: the parton densities in the virtual photon fi/ 7 *(x 7 * , P 2 , fx F „), 



for which up to the present there is only very little experimental information. This 
allows the study of the virtual photon structure by measuring jet cross sections. 

The contribution to the dijet cross section from resolved processes is expected 
to decrease relative to the contribution from the direct processes as P 2 — > fi 2 ; this 
means that the partonic content of the virtual photon becomes suppressed as P 2 
increases. Experimentally, resolved and direct processes are separated by using the 
variable x° BS , 

x obs = 1 (E ^i e -^«» + E^e-***), (1) 

where rf^ 1 ' 2 is the jet pseudorapidity and E e is the positron beam energy. x° BS 
measures the fraction of the virtual photon energy invested in the production of 
the dijet system and it is well defined at all orders in pQCD. For resolved (direct) 
processes, x° BS takes low (high) values. 

The measurements presented here were performed with the ZEUS detector at 
HERA. During 1995 to 1997 HERA operated with positrons of energy E e = 
27.5 GeV and protons of energy E p = 820 GeV. 



DIJET CROSS SECTIONS IN THE LABORATORY 

FRAME 

The x° BS distribution displays a high sensitivity to the virtual photon structure. 
Figure 2 shows the x° BS distribution for a sample of dijet events in three regions 
of P 2 . In all the P 2 regions studied, a significant resolved contribution is needed in 
the Monte Carlo to describe the data for x° BS < 0.75. Direct processes alone (the 
shaded area in figure 2) cannot describe the data [10]. 

Dijet cross sections as a function of x° BS have been measured for jets found using 
the longitudinally invariant kx cluster algorithm in the inclusive mode [11]. The 
measurements have been performed in the kinematic region given by 0.2 < y < 0.55 
and in three regions of P 2 : P 2 R3 (the photoproduction regime in which the 
scattered positron is lost in the beam pipe, figure 3a), 0.1 < P 2 < 0.55 GeV 2 (the 
intermediate P 2 region, in which the scattered positron is detected in the beam 
pipe calorimeter, figure 3b) and 1.5 < P 2 < 4.5 GeV 2 (the low P 2 DIS regime, 
in which the scattered positron is detected in the uranium-scintillator calorimeter, 
figure 3c). The events are required to have at least two jets with E T et > 5.5 GeV 
and —1.125 < i]^ et < 2.2. The data show that the shape of the measured cross 
section changes markedly with increasing P 2 : the cross section for low x^ BS values 
decreases faster than at high x° BS values. 

The measurements have been compared to several Monte Carlo models using 
various parametrisations of the photon parton densities (PDFs). The GRV [12] 
and WHIT2 [13] photon PDFs have been extracted for real photons and have 
no P 2 dependence. On the other hand, the SaS ID [4] photon PDFs consist of 
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FIGURE 2. The uncorrected x° BS distribution (black dots). Monte Carlo simulations using 
HERWIG are shown for comparison. 



two components with different P 2 dependence: the hadronic part decreases as 
~ ( m p/( m p + -P 2 )) 2 w ith increasing P 2 , while the P 2 dependence of the anomalous 
component goes like ~ log(/i 2 /P 2 ). The predictions of HERWIG [14] using the GRV 
PDFs do not describe the data, as expected since they have no P 2 dependence. The 
prediction based on WHIT2 provides a reasonable description of the data in the 
photoproduction and in the intermediate P 2 regimes. The predictions using SaS ID 
agree with the data at high x° BS in all regimes. The DIS Monte Carlo LEPTO [15] 
is compared to the data in the DIS regime. It underestimates the data at low x° BS , 
i.e. parton shower effects are not the sole contribution to this region of x° BS and 
P 2 . 

The low E 3 ^ 1 region is affected by contributions from the presence of a possible un- 
derlying event, which may mask the P 2 evolution of the virtual photon. Therefore, 
measurements of the dijet cross sections have been performed at higher transverse 
energies. For these measurements, two jets were required with E^ 1 ' 1 ^ > 7.5(6.5) 
GeV and —1.125 < rf et1 ' 2 < 1.875. The cross sections are shown in figure 4. At 
low x° the cross sections are also observed to decrease faster with P 2 than at 
high x® BS . The predictions from HERWIG based on the GRV PDFs are in good 
agreement with the data in the two lowest P 2 ranges, but fail at higher P 2 . On the 
other hand, the predictions based on SaS ID agree with the data at high P 2 but 
underestimate the data at low and intermediate P 2 at low x® BS . The predictions 
from LEPTO at high P 2 underestimate the data. 
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FIGURE 3. The measured dijet cross section da/dx^ BS (black dots) in the laboratory frame. 
The inner error bars represent the statistical errors of the data, and the outer error bars show 
the statistical errors and uncorrelated systematic uncertainties added in quadrature. The shaded 
band displays the uncertainty due to the absolute energy scale of the jets and the hatched band 
displays the uncertainty due to the modelling of the jet fragmentation. Monte Carlo calculations 
using HERWIG and LEPTO are shown for comparison. 



The P 2 evolution of the virtual photon structure was studied further by mea- 
suring the ratio of the dijet cross section for x° BS < 0.75 (enriched in resolved 
processes) to the dijet cross section for x° BS > 0.75 (enriched in direct processes). 
The ratio (see figure 5) falls steeply with P 2 which may be interpreted as the sup- 
pression of the resolved photon component as the virtuality of the photon increases. 
The prediction of HERWIG based on the GRV PDFs is constant, as expected from 
a photon structure without P 2 dependence. The prediction of SaS ID decreases 
with P 2 but lies below the data at low P 2 and the prediction from LEPTO shows 
that the contribution to the ratio from parton shower effects alone is not enough 
to explain the P 2 dependence of the data. QCD predictions have been calculated 
using the program JetViP [16] and different parametrisations of the photon PDFs. 
The predictions show sensitivity to the choice of PDF but lie well below the data. 
Hadronisation corrections (~ 20 — 30%) are insufficient to explain the discrepancy. 
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FIGURE 4. The measured dijet cross section da/dx^ BS (black dots) in the laboratory frame. 
Monte Carlo calculations using HERWIG and LEPTO are shown for comparison. 



DIJET CROSS SECTIONS IN THE j*p FRAME 

For P 2 ^> 0, the definition of x° BS (equation 1) is still valid. However, the 
photon remnant may have high transverse energy in the laboratory frame and be 
considered as a jet. The solution to avoid this problem is to transform to the 7*p 
frame. In such a frame the photon remnant has very low transverse energy and 
will not be mistaken for a jet emanating from the hard interaction. 

Figure 6a shows the measurements of the dijet cross section as a function of x® B 
in the 7*p frame [17]. The measurements have been performed in the kinematic 
region given by 0.2 < y < 0.55 and 0.1 < P 2 < 10 4 GeV 2 . The events are required 
to have at least two jets with E^j^ > 7.5(6.5) GeV and —3 < rfj?*' 2 < 0. Also for 
these measurements the shape of the measured cross section changes with increasing 
P 2 : the cross section for low x® BS values falls more rapidly with increasing P 2 than 
at high x° BS values. 

The predictions of HERWIG based on the SaS ID PDFs (see figure 6a) give 
a good description of the data in the high P 2 region but fail in the intermediate 
P 2 region. A resolved photon component is needed to describe the data up to 
P 2 ~ 49 GeV 2 . Above this value, the HERWIG prediction is dominated by the 
direct component and describes the data well. Therefore, for P 2 > (E^) 2 , direct 
processes alone are able to describe the data. 

The ratio of the cross section for x^ BS < 0.75 and x° BS > 0.75 also falls steeply 
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FIGURE 5. The ratio of dijet cross sections, a(x° BS < 0.75)/a(x° BS > 0.75), as a function of 
P 2 (black dots). Monte Carlo (HERWIG and LEPTO) and QCD (JetVip) calculations are shown 
for comparison. 



in the 7*p frame with increasing P 2 (figure 6b). The prediction of SaS ID shows a 
decrease with P 2 but lies below the data in the whole P 2 range. The ratio of the 
data to the SaS ID prediction (inset in figure 6b) has a constant value of ~ 1.3. 
This indicates that the resolved component suppression included in these PDFs 
is in agreement with the data, but they underestimate the fraction of resolved 
component by ~ 30%. 



SUMMARY AND CONCLUSIONS 

Dijet cross sections as a function of x® BS have been measured in the laboratory 
and 7*p frames in the kinematic region given by 0.2 < y < 0.55 and < P 2 < 10 4 
GeV 2 for jets found using the longitudinally invariant fcr cluster algorithm in the 
inclusive mode. The x° BS dependence of the measured cross sections changes with 
increasing P 2 : the cross section for low x° BS values decreases faster than at high 
x obs va i ues Tlig predictions of HERWIG based on the SaS ID photon PDFs 
describe the data for 1.5 < P 2 < 10 4 GeV 2 but fail in the region 0.1 < P 2 < 0.55 
GeV 2 . A resolved photon component is needed to describe the data up to P 2 ~ 49 
GeV 2 , i.e. where the virtual photon is probed at a scale comparable to the hard 
interaction scale (E^) 2 ~ (7 GeV) 2 . The ratio of the dijet cross section a(x° BS < 
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FIGURE 6. (a) The measured dijet da/dx° BS (black dots) in the j*P frame, (b) The ratio of 
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dijet cross sections, a{x^ uti < 0.75) /<r(x° BS > 0.75), as a function of P 2 (black dots). Monte 
Carlo calculations using HERWIG are shown for comparison. 



0.75)/a(x° BS > 0.75) decreases with P 2 ; the predicted P 2 dependence of the ratio 
agrees with the data, but it underestimates the fraction of the resolved component 
by ~ 30%. This result can be interpreted in terms of a resolved photon component 
that is suppressed as the photon virtuality increases. 
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